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X-ray tomography can provide structural information of whole cells in close to their native state. 
Radiation-induced damage, however, imposes a practical limit to image resolution, and as such, a choice 
between damage, image contrast, and image resolution must be made. New coherent diffractive imaging 
techniques, such Fresnel Coherent Diffractive Imaging (FCDI), allows quantitative phase information with 
exceptional dose efficiency, high contrast, and nano-scale resolution. Here we present three-dimensional 
quantitative images of a whole eukaryotic cell by FCDI at a spatial resolution below 70 nm with sufficient 
phase contrast to distinguish major cellular components. From our data, we estimate that the minimum 
dose required for a similar resolution is close to that predicted by the Rose criterion, considerably below 
accepted estimates of the maximum dose a frozen-hydrated cell can tolerate. Based on the dose efficiency, 
contrast, and resolution achieved, we expect this technique will find immediate applications in tomographic 
cellular characterisation. 

X-ray Coherent Diffractive Imaging (CDI) is currently under intensive development for a wide range of 
applications in the materials 1 " 3 and biological sciences 4 ' 5 , due in part to its ability to provide spatial 
resolution beyond the limits imposed by image-forming optics 6 . The short wavelength and element specific 
interaction of X-rays with matter may be exploited to characterise intracellular features 7 " 10 with minimal sample 
preparation. In applications involving CDI of biological samples, the ultimate resolution limit will generally be 
determined by radiation damage rather than the coherent power of available light sources 1112 . In the context of 
cellular imaging, it is therefore crucial to optimise the sensitivity and dose efficiency of the technique. 

Historically, CDI could only be applied to small isolated samples due to the limits imposed by sampling of the 
diffracted wavefield 1314 . Recent modifications to the standard CDI algorithm including Fresnel CDI (FCDI) 1516 
and ptychography 17 overcome the issue, allowing the CDI method to be applied to objects of arbitrary size. FCDI 
utilises illumination with known phase curvature to provide additional redundancy, providing more reliable, 
unique image reconstruction and increased tolerance to partial coherence from a single probe position 151819 , 
while ptychography involves collecting data from multiple overlapping probe positions, significantly increasing 
image quality 20 . Ptychographic plane-wave CDI has been applied to biological samples 20 " 22 and has demonstrated 
a very good compromise between dose (2.5 X 10 8 Gy) and resolution (60 nm) when applied to three-dimensional 
imaging of a bacteria cell 23 . However, the required imaging dose predicted by the Rose criterion for a resolution of 
60 nm is only 10 6 Gy 11 . Ptychographic FCDI has also been demonstrated for several applications and has been 
shown to deliver the advantages of both of the techniques it combines without significantly adding to the 
experimental complexity 24 " 26 . In particular, this method has been demonstrated with known test samples in 
two- dimensions 24 and three-dimensions 25 . 

To demonstrate the viability of ptychographic FCDI for whole-cell three-dimensional imaging, a red blood cell 
(RBC) infected with the malaria parasite Plasmodium falciparum at the trophozoite life-cycle stage was used as a 
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Figure 1 | Three-dimensional rendering of the reconstructed phase of the 
sample shown at two orthogonal angles. The whole infected red blood cell 
(blue), with the parasite (red) and the parasite's digestive vacuole (black) 
are illustrated in orthogonal views a) and b). Views c) and d) show the 
surface renderings of the parasite and the digestive vacuole corresponding 
to a) and b) respectively to highlight the parasite's exomembrane 
structures. These are areas where the parasite extends into the host cell and 
likely correspond to components of the exomembrane system such as the 
tubulovesicular network. The scale bar is equal to 1 urn. 



sample. During the blood stage of the P. falciparum life cycle, the 
parasite catabolizes hemoglobin breaking it down to amino acids 
and releasing free heme, which is toxic to the organism. The malaria 
parasite detoxifies the heme by crystallizing it to form malaria pig- 
ment or hemozoin 27 . The variation in structure and density between 
the high-density hemozoin and the low-density hemoglobin pro- 
vides an ideal test for high contrast biological imaging in three 
dimensions. From the data presented here, the lower limit of deliv- 
ered dose required for three-dimensional images at a defined con- 
trast and the achieved spatial resolution of 70 nm is evaluated. 
Additionally, by recovering a three-dimensional map of the com- 
plex transmission function of the sample, the complex refractive 
index in each volume element of the cell may be determined, which 
facilitates the identification of ultracellular features through image 
segmentation. 



Results 

Ptychographic FCDI diffraction data were collected from a total of 79 
individual projections spanning 124.5 degrees at an X-ray energy of 
2.5 keV. At each projection angle, the sample was scanned through 
the beam to collect data from nine overlapping areas of the sample. 
Supplementary Figure SI depicts an example set of diffraction data. 
The total accumulated imaging dose is less than 2.35 X 10 8 Gy. 
Retrieval of the complex sample transmission function for each pro- 
jection was performed using the NADIA software library (NADIA 
Software Project http://www.coecxs.org/joomla/index.php/research- 
and-projects/nadia-software-project.html) and the 3-D volume 
reconstruction was carried out by algebraic reconstruction tomo- 
graphy (ART) based on the Kaczmarz algorithm 28 . Segmentation 
was performed based on the measured refractive index parameters 
described below and volume rendering was done using the Amira 
software package. Three-dimensional images of the phase recon- 
struction of the transmission function of the sample at two ortho- 
gonal angles are shown in Figure 1 a) and b). 

Three regions are indentified corresponding to the host RBC (ren- 
dered in blue), the intraerythrocytic parasite (rendered in red), and 
the hemozoin- containing digestive vacuole (rendered in black). 
The advantages of the 3D volume reconstruction process are evi- 
dent by comparing projections of the reconstructed 3D volume 
(Figure 1) with an individual two-dimensional ptychographic image 
(Figure 2b). The features are best appreciated by examining a rota- 
tion of the rendered 3D model (Supplementary Video). The parasite 
has an amorphous shape as described previously for trophozoite 
stage parasites 29 . Regions of the surface of the parasite show exten- 
sions into the RBC cytoplasm, highlighted in the surface rendering of 
the parasite and digestive vacuole with the host RBC removed in 
Figure 1 c) and d). These likely correspond to components of the 
exomembrane system of the parasite, such as the tubulovesicular 
network 30 . The host RBC is deformed from the classical biconcave 
disc shape of uninfected RBCs reflecting the substantial reorganisa- 
tion of the host cell by the intraerythrocytic parasite 31 . 

Quantitative measurements of the complex refractive index taken 
from the tomograms were used to isolate the three different compart- 
ments comprising the host red blood cell, the parasite, and the para- 
site's digestive vacuole. Values for the decrement from unity of the 
real part of the complex refractive index, 5, were measured to be 2.2 
± 0.7 X 10" 5 , 3.9 ± 0.4 X 10" 5 , and 5.2 ± 0.4 X 10" 5 respectively. At 
2.5 keV, these values are in agreement with the tabulated values 32 
for hemozoin, a significant proportion of the digestive vacuole 
(C 3 4H3 2 0 4 N4Fe 1 , density 1.45 g/cm 3 ) 33 with a value of 5 = 5.12 X 
10" 5 , and haemoglobin, contained in the red blood cell 
(C 738 H 1166 N 812 O203S2Fe 1 , density = 0.38 g/cm 3 ) 34 with a value of 5 




Figure 2 | Line profiles through one projection angle with exposure times of 20, 10, 5, 2, and 1 second per ptychographic point (a). The results indicate 
little change for between 20 and 5 seconds, with 1 and 2 seconds of exposure per point (dotted black and dashed green lines respectively), leading to 
significant variations in the reconstructed phase values (see inset). From this result we can extrapolate that 5 seconds of exposure time per point would 
yield images of similar quality, with one fifth of the present acquisition time using the experimental configuration described in the methods section. 
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Figure 3 | Resolution of the reconstructed volume: The PTRF 
(Ireconstructed/Imeasured) as a function of the spatial frequency. Also 
displayed are the theoretical MTF with 70% efficiency for Abbe resolutions 
equal to 55 and 70 nm. The measured PTRF falls between these values, 
indicating resolution of better than 70 nm in the image. The reconstructed 
image was taken as a maximum intensity projection through the 3D 
volume and was compared to diffraction data from the same rotation 
angle. At low spatial frequencies, the beamstop on the zone plate masks the 
diffraction data and was therefore not included. 

= 1. 35 X 10~ 5 . The parasite itself contains many protein based 
compounds, and was therefore modelled as empirical protein 
(CaoHsoOxoSx), however, the density of the parasite changes during 
the life cycle stage 34 . For the given composition and life cycle stage, 
the measured value of 5 for the parasite implies that the density is 
equal to 1.1 g/cm 3 . The higher density of the cytoplasm of this tro- 
phozoite stage parasite compared with the host cell may reflect the 
high density of ribosomes at this stage 35 . 

While the numerical aperture of the diffracted signal (related to 
the scattering strength of the sample 11 ) sets a theoretical upper limit 
on the resolution, it is the extent to which this data can be accurately 
phased which practically determines the image resolution in CD I 8 . 
One method that is often used in determining the spatial resolution 
of the reconstructed images is to calculate the Phase Retrieval 
Transfer Function (PRTF). The PRTF is estimated by calculating 
the ratio of the intensity of the phased diffraction pattern to the 
measured intensity 2 . In this case, a projection through the transmis- 
sion function of the entire reconstructed volume was created, which 
was propagated to the detector plane using the Fresnel propagator 36 . 
This was compared with the measured intensity data for the same 
projection angle after the incident illumination had been subtracted. 
This measurement of the resolution of the entire reconstructed 
volume therefore takes into account slice artefacts and artefacts aris- 
ing from the limited angular scan range. As a comparison, the PRTF 
is shown alongside the theoretical modulation transfer function 
(MTF) for comparable conventional optics with 70% efficiency 
and an Abbe resolution equal to 55 and 70 nm in Figure 3. We 
conclude from this analysis that the spatial resolution in the three- 
dimensional image lies between 55 and 70 nm. 

Putkunz et al 26 showed that by using many overlapping probe 
positions in FCDI, more accurate phasing of the diffraction data is 
possible while the total dose delivered to the sample can be greatly 
reduced without reducing the quality of the reconstructed image. To 
examine this result in the context of cellular tomography, we recon- 
structed images from diffraction data obtained with exposure times 
ranging from 1 s to 20 s per position in a ptychographic series at one 
particular projection angle. Line profiles through the reconstructed 
phase are shown in Figure 2, together with an example reconstruc- 
tion in which their position in the sample is indicated. The similar- 
ities in the profiles show that the quality of the image is consistent 
between the datasets with exposure times of 5 s, 10 s, and 20 s, with 
an effective dose of between 6.0 X 10 5 and 2.4 X 10 6 Gy, in line with 
previously reported results 26 . This provides confirmation that 3D 
cellular diffractive imaging can be realised using ptychographic 



FCDI with a lower dose than (full-field) transmission X-ray micro- 
scopy (TXM) due to the inherent phase sensitivity of the former and 
the lack of attenuating optics between the sample and detector 37 . We 
note that Scanning Transmission X-ray Microscopy (STXM) also has 
no post-sample optics. High resolution STXM tomography is tech- 
nically challenging, quantification generally demands increased 
dose, and the small depth of focus of zone plate optics limits sample 
thickness 39 . FCDI is also limited by sample thickness, albeit though 
different mechanisms, but this can be partially overcome by tailoring 
the experimental geometry to the sample thickness. Additionally, 
unlike STXM, the resolution of CDI is not limited by the NA of 
the optics 15 . Based on the analysis above, we can extend the hypo- 
thetical dose reduction to the whole dataset, yielding a 3D recon- 
struction with similar quality with a dose of 4.7 X 10 7 Gy, 
approaching the minimum required dose for imaging at the esti- 
mated spatial resolution reported by Howells et al. of approximately 
10 6 Gy 11 . This result compares favourably in terms of dose and reso- 
lution to transmission X-ray microscopy (TXM), where a dose of up 
to two orders of magnitude greater may be required to achieve a 
similar resolution 38 . Our analysis shows that sub -70 nm imaging 
resolution can be achieved for thick, whole cells by ptychographic 
FCDI with an X-ray dose considerably below the tolerable maximum 
determined for frozen hydrated specimens 11 . While image resolution 
for such samples is expected to be ultimately limited to about 10 nm 
by radiation induced damage 11 , it is clear that coherent X-ray flux 
available from synchrotron based light sources will impose a prac- 
tical limit on the achievable resolution 12 . Recent advances how- 
ever, which include the coherent properties of the source 19,39 
combined with emerging detector technologies 40 , offer the possibility 
of allowing 3D cellular FCDI data to be acquired in less than one 
minute. 

Discussion 

We have demonstrated whole-cell ptychographic FCDI tomography 
with high contrast and dose efficiency to sub 70 nm resolution. High 
contrast in the three-dimensional map of refractive index made it 
possible to identify structures internal to the cell through segmenta- 
tion, such as the ribosome- and membrane-rich parasite volume with 
a higher density than the host cell, but a lower density than the 
parasite's digestive vacuole (which houses crystals of hemozoin). 
Additionally, the analysis shows that a similar quality reconstruction 
can be achieved with a dose approaching the minimum dose required 
for imaging at the resolution achieved in the present work 11 , two 
orders of magnitude less than TXM 38 , and leading not only to lower 
damage sustained by the sample, but also decreasing the acquisition 
time substantially. Combined with recent advances in partial coher- 
ence and pink beam methods 39 these results offer the possibility to 
reduce the acquisition time to several seconds, allowing the possibil- 
ity of statistical numbers of biological samples to be imaged in a 
reasonable time frame. This therefore increases the value of the tech- 
nique to the biological sciences. Extending to frozen hydrated sam- 
ples in the water window, a similar result with respect to resolution 
and dose should be achievable 11 . Indeed, the increased available con- 
trast between carbon containing proteins and the background 38 
could result in a lower dose while maintaining high contrast images. 

Methods 

Our experiment was conducted at Beamline 2-ID-B at the Advanced Photon Source 41 
using an in vacuo imaging endstation 42 that permits a sample to be rotated and 
translated through a diverging beam formed via a Fresnel zone plate (FZP). We used a 
2.5 keV X-ray beam to coherently illuminate a 160 um FZP with nominal outer zone 
width of 30 nm, and placed the sample within the diverging beam. At this energy the 
focal length is 9.81 mm. All but the first order diffracted beam was blocked by a 
combination of a 40 um central beam stop at the zone plate and a 10 um order sorting 
aperture (OS A) positioned near the focal plane. 

P. falciparum was cultured in vitro and mature stage-infected red blood cells were 
harvested and fixed in 0.01% paraformaldehyde before being freeze-dried for 30 
minutes. The sample was mounted on a gold electron microscopy grid covered with 



SCIENTIFIC REPORTS | 3 : 2288 | DOI: 1 0.1 038/srep02288 



3 



formvar. While this assisted the sample location process, the tilt range was limited to 
124.5 degrees. The sample was initially placed as close as practical to the focus of the 
diverging beam to maximize the number of incident photons. At high tilt angles, 
geometrical constraints led to the sample being positioned further downstream. Over 
the entire tilt range, the sample was placed between 375 um and 650 urn downstream 
of the focus where the illumination diameter was between 6 urn and 1 1 urn. This 
variation has the benefit of introducing a form of phase- diversity into the three- 
dimensional reconstruction, which has been shown to reduce artefacts due to the 
structure in the illumination associated with the zone plate 26 . Ptychographic dif- 
fraction data in a 3 X 3 array with a 1.5 um step size were collected at 84 projection 
angles uniformly spaced over 124.5 degrees. At each ptychographic position, 10 one- 
second exposures were summed to increase the signal to noise ratio and dynamic 
range in the data, while every third projection angle consisted of the sum of 50 one- 
second exposures per ptychographic position, which assists the 'bootstrapping' 
process 43 . An additional 15 projections were taken at the greatest incident angles to 
reduce tilt increment to ameliorate the effect of the projected increase in substrate 
thickness. Five additional projections were taken about normal incidence. Of the 104 
projections measured, 25 were discarded due to poor reconstruction quality, typically 
due to uncontrolled sample movement during data acquisition, resulting in the total 
accumulated imaging dose of less than 2.35 X 10 8 Gy. 

Diffraction data were collected on a cooled charge-coupled device detector 
(Princeton MT-MTE) cropped to 1 150 X 1 150 square pixels each with a 13.5 um side 
length placed 83 cm downstream of the sample. The corresponding pixel size in the 
sample plane due to the numerical aperture of the detector is 26 nm. The actual 
resolution of reconstructed image depends on the signal-to-noise ratio at large angles 
in the diffraction pattern. Based on Abbe theory 44 and the numerical aperture (NA) of 
the detector (defined here to be one half of the detector acceptance angle) we find a 
theoretical resolution limit of T = 0.82^/NA ~ 42 nm. 

Each projected angle is reconstructed individually enforcing phase and amplitude 
complex constraints 45 using the NADIA software library (NADIA Software Project 
http://www.coecxs.org/joomla/index.php/research-and-projects/nadia-software- 
project.html). To increase the speed of convergence of the reconstruction algorithms, 
the transmission function from each projection was used to bootstrap the subsequent 
projections, leading to a faster, more reliable convergence for each projection 43 , as 
illustrated by Peterson et al. 25 . Using this method, projections with 50 exposures per 
ptychographic position were reconstructed with 2000 iterations of the ER phase 
retrieval algorithm, while the remaining projections were reconstructed with 400 
iterations. The individual projections were coarsely aligned before sinograms were 
created. Jitter in the sinograms was corrected using the iterative method that does not 
require fiducial markers inserted within the specimen described by Abbey et al. 
(Abbey, B. et al. Iterative Marker- Free Alignment Method for Low Quality 
Tomography Data, unpublished data). The tomograms were then reconstructed 
using the algebraic reconstruction tomography (ART) algorithm based on the 
Kaczmarz algorithm 28 . Due to the relatively high number of projections, qualitative 
convergence of the algorithm was achieved within 10 iterations. Segmentation based 
on the measured refractive indices of each reconstructed voxel, as described in the 
results section, and volume rendering was done using the Amira software package. 
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